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Abstract: In this letter, I solve a model for the dynamics of vortices in a decaying 
two-dimensional turbulent fluid. The model describes their effective diffusion, and the 
merging of pairs of vortices of same vorticity sign, when they get too close. The merging 
process is characterized by the conservation of energy and of the quantity Nr n , where 
r is the mean vortex radius, and iV their number, n = 4 corresponds to a constant 
peak vorticity, and n = 2 to a constant kurtosis. I found the scaling laws for various 
physical quantities (r, enstrophy, kurtosis...), and for instance, it is shown that N ~ 
(Vln(t)) 3«-4 for n > 2, and N ~ t~ L for n = 2, in good agreement with extensive 
numerical simulations. I also discuss some recent experiments in view of these results. 
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Recent experimental [1-2] and theoretical [3-7] works have emphasized the importance 
of coherent vortex dynamics during the decay of turbulence in d = 2. This process 
consists mainly in three stages: during an initial transient period, the fluid self-organized 
and a network of coherent vortices appears. Once the coherent vortices have emerged, 
vortices disappear essentially through merging of same sign vortices, such that their 
number decreases and their average size increases, in a process somewhat reminiscent of 
a coarsening stage [8] [9]. Finally, when only one (or very few) dipole is left, it decays 
diffusively, due to the finite viscosity. 

In order to describe this second stage of the dynamics of decaying turbulence, the 
authors of [6] [7] introduced the following simple model. The starting point lies in the well 
known fact that the evolution dynamics of a set of far apart vortices is a conservative 
process, such that vortex centers positions {Xj} evolve according to Kirchoff laws [10], 

r/X- 

r l ^ = v ! xw, w = -^r i r i in|x i -x i | (i) 

where the circulation I\ = nrjujj is the area of the i-th vortex multiplied by its average 
core vorticity. In addition to the dynamics, one must describe the dissipation process 
when some vortices get too close. Experiments [1] [2] or numerical simulations [5-7] have 
consistently shown that this happens mainly through a complicated merging process 
between vortices of same vorticity sign, while dipoles form a very stable state [2]. In 
[6-7], this was modelized by assuming that when two vortices with radii r\ and r<i are at 
a distance less than a dj^r^r^), they instantly merge. The authors of [6] found that a 
consistent form for this minimal approach distance is d c (r\, T2) ~ 2.592r2 + 0.609r^/r2, 
for r\ < T2- In the present paper, the claim is made that the precise form of d c is not 
important as far as the critical properties of the model are concerned, provided d c is of 
order of the mean vortex radius r. When two vortices of same sign merge, we still need to 
define the properties of the resulting vortex. Motivated by their numerical results for the 
full Navier-Stokes equation, the authors of [5-6] were lead to impose local conservation 
of energy and a time independent typical peak vorticity u>. This was also found to be 
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consistent with the experiment by Tabeling et al. [1]. Since the energy of a vortex is 
of order T 2 ~ c<j 2 r 4 , the new vortex has a radius r% satisfying, r| = r 4 + r|. On the 
other hand, following the observation of the experiment described in [2] (see below), 
we could alternatively impose that the energy and the total area occupied by vortices 
remains constant, which is essentially equivalent to have a constant kurtosis K. The new 
merging rule is then, r\ = r\-\-r\, and u must decay in order to keep the energy constant. 
More generally, we will treat the general case of a dynamics with the conservation of the 
local "charge" qi = rf. 

We now derive the scaling relations between the number of vortices N, their mean 
radius r, their typical distance /, the ratio Z/E (Z is the enstrophy), the kurtosis K, and 
uj/y/E. For large time, we expect that r ~ t v [1-2] [5-7] [11] [12]. The assumed conservation 
law, and the fact that E ~ iVcu 2 r 4 , Z ~ Nu 2 r 2 , K ~ iVr 2 , leads to, 

-~t- 2 -, K~ t v- n >, 4L~H 2 -§)- (2) 

E ' VE 

which was obtained in [5-6] for the case n = 4. We need n > 2 to keep / > r, while 
for n < 2, one can show that iV ~ iVo exp(— t/to) [9]. We see that the problem now 
amounts to determine the unknown exponent v which should depend explicitly on the 
conservation law considered in the dynamics. Ultimately, one should be able to infer 
on physical grounds (role of initial conditions [7], viscosity, or d = 3 effects [2]) which 
effective conservation law, or n, is relevant in experiments or numerical simulations of 
(quasi-) two-dimensional decaying turbulence. I now focus on the relation between v 
and n, and sketch a rather simple and physical derivation for it, while a more elaborate 
solution is possible (see conclusion). This solution is based on the mapping of the original 
problem on a solvable dynamics for a self-organized system [13-14]. This, maybe, gives 
a more solid justification for the description of the coherent vortex structure in terms of 
a self-organized system [15], as mentioned by some authors [6]. 

The first step of the solution uses the claim that the motion of vortices in the absence 
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of merging is effectively diffusive, which should result from the chaotic nature of their 
motion. From Kirchoff equations, and using the fact that the total circulation is zero, the 
mean square velocity of vortex centers is, v 2 = -j^r^ ~ Noj^r 4 " J^ 1 R ^ ~ Eln(L/r), 
where the upper cut-off is the linear size of the system L, and the lower one is the mean 
vortex radius r. Note that v 2 is strictly constant with the correct definition of energy 
E ~ iVu; 2 r 4 ln(L/r), but decays slowly with the definition E ~ iVu; 2 r 4 considered in the 
simulations of [6] and also here, for the moment. The diffusion constant of a vortex is 
then of order D ~ v /u, which varies in time if u is not constant. We can now consider 
separately the two gas of vortices with positive and negative vorticity, since the effect 
of one on the other is taken into account in the diffusion constant, and since mergings 
only involve same sign vortices [2] [5-7] . It now remains to treat correctly the apparently 
complicated merging condition. This can be done by mapping the problem to an effective 
dynamics on a lattice of suitable constant a, and by choosing consistently a time step 
t = a? /D. If one takes a = r, the merging process is correctly described by saying that 
two charges add together when they hop on the same lattice site. Indeed, in this case their 
mutual distance is then less than the mean vortex radius, which is in accordance with the 
kind of merging conditions introduced above. We have thus separated the dynamics in 
two processes: (i) free diffusion of far apart charges (vortices), with a diffusion constant 
D ~ f 2 /cu, and (ii) addition (merging) of charges hopping on the same site (vortices 
at a distance less that a = r). What makes the problem solvable is the fact that the 
dynamical equations for the time evolution of the q^s are now linear: 

Qi(t + t)= J2 w j^iQj(t), w j^i = or 1, J2 w j^i = 1 ( 3 ) 

\i-j\=a i 

where the Wj^s are random variables, and the last condition expresses that the vortex 
and its associated charge on site j hops on one of its neighboring sites (with equal 
probabilities). This problem can be solved exactly [13-14], even in the case when one adds 
a random charge 7j(£) to the left hand side of Eq. (3), which could simulate the driving 
noise in fully developed turbulence [9]. For the purpose of this letter, we do not need the 
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full solution of Eq. (3) , but only need to know that the upper critical dimension beyond 
which mean-field theory is exact is d = 2 [16] [14]. Assuming no correlation (mean-field), 
the number of collisions at time t is proportional to iV 2 . Then, ^ ~ —(a/L) d N 2 /r : 
which immediately leads to N ~ (L / a) d (t / r) -1 , where L is the linear size of the system 
and d > 2 is the space dimension. For d = 2 and for n > 2 (l/r then diverges when 
t — > oo and vortices are point-like objects), one can show exactly [17] the existence of 
logarithmic corrections in this kind of problems, such that N(t) ~ L 2 \n(t / t) / (Dt) . 

In order to illustrate the quantitative effect of such a correction on the numerical 
determination of the decay exponent, I show in fig. 1 numerical simulations of Eq. (3) 
with a constant D = a? jr on a 1500 x 1500 square lattice with all sites initially occupied, 
and up to t = 10 5 r. The best fit by the law N(t) = N (l + t/t ) a (as used in [6]) gives 
a remarkably linear fit with a = 0.90 ± 0.01 even after five decades in time (see fig. 1-a)! 
Of course, the plot of tN(t) vs ln(£) immediately reveals the logarithmic correction (fig. 
1-b). This remark motivates my next comment: there have been some claims that the 
vortex dynamics in decaying turbulence presents some similarities with the coarsening 
dynamics of the XY model [18], for which the number of spin vortices was apparently 
seen numerically to decay as iV ~ t~ - 75 [18-19], as in the Kirchoff vortex model [6]. 
Recent more extensive simulations [20], have shown that this apparent behavior is due 
to an important logarithmic correction, and that iV ~ ln(f)/£, as expected from the 
theory of coarsening systems with non conserved order parameter [8] [20]. Since the 
anomalous exponent £ seen in the XY model seems to be an artifact of unanticipated 
logarithmic corrections, it is natural to challenge the observation of such exponents in 
the field of decaying turbulence. Notice however that the vortices in the XY model 
behave in a qualitatively quite different fashion as in a turbulent fluid [9]. First of all, 
they are quantized, and after a short transient time, only vortices with winding number 
±1 remain [20] [8]. Contrary to the merging dynamics of fluid vortices already described, 
that of the XY model only involves the annihilation of opposite (quantized) vortices. 
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After these preliminary remarks, I end the derivation of the scaling exponent u, and 
confirm the role of logarithmic corrections. Using the scaling equations Eq. (2), and the 
above result for N(t), one easily finds, r ~ (v^t/ ln(t))^, with v = p^z^y, since v does 
not behaves faster than a logarithm in time. 

For the first case of a dynamics with conserved energy and constant ui (n = 4) 
studied in [5-6] [7], the result is r ~ (v 2 t/ hit) 1 / 4 , with v = \n(t en d/t), and where t en( i 
is of the order of the time for which only one (pair of positive and negative) vortex is 
left in the sample. We thus find v = 1/4, or in the language of [5-6], £ = &v = 1. 
The authors in [6] found £ = 0.72 ± 0.02. I will now show that the discrepancy with 
the present analytical result is certainly due to the two logarithmic corrections, the 
first one in v, and the second one coming from the critical d = 2 corrections, which 
both tend to lower the effective exponent £. In fact, it would have been more correct 
to consider the definition E ~ iVu; 2 r 4 ln(L/r) [6], taking into account the long range 
interactions between vortices, leading to the local conservation of r 4 ln(L/r). This implies 
that the typical vortex square velocity is now constant. Still, there remains the most 
important logarithmic correction (due to d = 2), which in principle, should be present 
in the numerical solution of the full Navier-Stokes equation, for which the authors of [6] 
also found £ ~ 0.70 ~ 0.75. In fig. 1, we show the result of numerical simulations of the 
problem of diffusing vortices merging when their distance is less than d c (ri, 7*2) = r\ + r2, 
and with local conservation of r 4 , but with a constant v 2 (and D since u is constant for 
n = 4), in order to avoid the introduction of extra logarithmic corrections (present in [6]). 
One finds £ = 4u = 0.84 ± 0.01 (fig. 1-a), from a fit to the form N(t) ~ N (l + t/t )~^ 
as used in [5-6], even when using much bigger samples than in [6-7] (Nq = 20000 instead 
of iVo ~ 400 in [6] and Nq = 300 in [7]), which is allowed by the simpler nature of the 
model. However, for this model, the exact solution presented here predicts N(t) ~ \n(t)/t, 
which is illustrateded on fig. 1-b by plotting tN(t) vs ln(£). I now present numerical 
simulations of the full model introduced in [6-7], where the vortices evolve according to 
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Kirchoff dynamics (Eq. (1)), and for which the same merging criterion as in [6] as been 
used. The calculations are much heavier than for the simple diffusing vortex model, which 
drastically limits the number of vortices and integration time that one can simulate. In 
order to obtain the properties of the scaling regime, the authors of [6] introduced a 
clever renormalization scheme: they have simulated the decay of the system from 400 
to 100 vortices, and used the final configuration to generate a new initial condition with 
400 vortices, by copying it four times. They repeated this renormalization scheme until 
they obtained an apparent scaling regime and enough statistics. The drawback of this 
method is that it only allows the precise measurement of the decay of the system during 
less than a decade in time, during which the number of vortices is divided by a factor 
4 (in fact less [6] [9]; see also fig. 1). They did not find a pure power law decay, but 
found a good fit to N(t) ~ JV (1 + t/t )~^, with f « 0.72. This form simulates the 
effective (sharp) increase with time of the exponent £, as they found quite large values 
for to ~ t en( }/3 (t en d ^ s the final time of the simulation). This scaling was shown in [11] 
to be consistent with the assumption that the population of vortices is at all time in 
local Boltzmann equilibrium. Note however that this assumption would be wrong in a 
coarsening dynamics (for instance for the XY model), for which the correlation functions 
during the dynamics are unrelated to the equilibrium ones [8] [9] [20] . Such a scaling with 
roughly the same effective n = 4 and £ as 0.7 was also observed experimentally in [1]. 
As far as the direct numerical simulation of the Navier-Stokes equation is concerned, the 
authors of [6] claim to see the same behavior as described above, while some other groups 
discuss the existence of such growth exponents [21], or their values [22]. The simulations 
presented in fig. 1 use the same method as in [6] , but with an initial and final number of 
vortices respectively equal to 1800 and 200 (instead of 400 and 100). The details will be 
given in [9]. The apparent exponent as extracted by the method of [6] is £ = 0.76 ± 0.02 
(fig. 1-a), systematically larger than in [6]. In [7], where smaller samples (N = 300) and 
shorter times than in [6] where considered, the authors found an even lower £ ~ 0.6. In 
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fig. 1-b, I show the plot of tN(t) vs ln(rj), which is found to be reasonably linear. Note 
that logarithmic corrections as found here due to the merging (coarsening) dynamics of 
point-like objects (r <C /) in d = 2 could also be generated by taking into account [9] the 
long range interaction between vortices as in the XY model [20]. Finally, the mean-field 
vortex radii distribution can be computed [9] using Eq. (3) and the results of [14] . Once 
radii are expressed in unit of the average radius r, the result is P(x) = 4Ax 3 exp(— Ax 4 ), 
where A = T(5/4). The distribution is shown on fig. 2. The agreement with numerical 
simulations is satisfactory, the slight discrepancy being probably due to the difference 
between the merging distance used in the simulations (d c = r\ + 7*2) and in the exact 
calculation (d c = r), and the fact that some small correlation (non mean-field) effects 
should be present in d = 2. The distribution for the Kirchoff dynamics is slightly wider 
[6], but d c is probably too different from the merging condition in the solvable model. 

We now study the case of the coherent vortex dynamics with conserved kurtosis. In 
other words, the area covered by the vortex cores remains constant. This phenomeno- 
logically reproduces the observed feature of the experiment presented in [2]. If we also 
impose conservation of energy, the maximal vorticity must then decay as u> ~ r _1 accord- 
ing to the relation E ~ iVr 4 u; 2 ~ r 2 c<j 2 , for constant K. Noting that D ~ v^/uj ~ r ~ / 
increases with time, one finds r ~ v^t, such that v = 1, with no logarithmic corrections 
(for constant v) since r ~ /, and thus vortices are no more point-like objects [9]. This 
exactly reproduces Batchelor phenomenology [12], in the framework of a microscopic 
theory. I have simulated Kirchoff's dynamics for n = 2, with initially N = 1600 vortices. 
After one renormalization step, iV is left decaying. In this case the dynamics is much 
faster than for n = 4, which permits to observe the decay up to iV ~ 0(1) (fig. 3). 
The fit to the form N(t) = Nq(1 + t/to) ~ 2 is excellent (fig. 3-a). Similar results are 
obtained (on more decades) for the diffusing vortex model with D ~ r, confirming v = 1 
(fig. 3). In the experiment described in [2], the obtained scaling relations are exactly 
of the form given in Eq. (2) with n = 2, but the experimental estimate for v (in the 
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first time decade) lies in the 0.18 ~ 0.26. In fact, a naive determination of 

v in a Kirchoff vortex simulation with initially Nq — 100 vortices, in the time range for 
which 0.2iVo < N < 0.8iV (as in [2]), leads to 2u = 0.54 ± 0.04 in surprisingly good 
agreement with experiment (see fig. 3 with Nq = 1600). Note that three-dimensional 
effects, leading to the diffusion of vorticity in the third dimension, could be at the origin 
of the constant area occupied by vortices in this experiment [2]. 

In conclusion, I have solved a simple model of diffusing vortices, with the effective 
diffusion constant derived from the Kirchoff dynamics and conservation laws. I have 
found long transient regimes (due to logarithmic corrections for n > 2) which, I argue, 
are at the origin of the smaller exponents observed in the literature. For a constant 
peak vorticity u (n = 4), vortices are effectively diffusive (D ~ v 2 /uj is then constant), 
while for a constant kurtosis K (n = 2), their motion is effectively ballistic (D ~ t). In 
the presence of a finite viscosity (diffusive dissipation), one could expect that the first 
dynamics is more stable than the Batchelor regime, the motion of ballistic vortices being 
more easily damped than that of diffusing ones. Note finally that the problem of diffusing 
vortices can be solved in many other situations: finite driving noise, initial long range 
distribution of vortex radii (as numerically studied in [7]),... This, and more details on 
the numerics, as well as a new field theory for d = 2 decaying turbulence in the spirit of 
the vortex model will be presented in [9] . 
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FIGURE CAPTIONS 



Fig. 1: Number of remaining vortices N(t) after a time t for the different models introduced 
in the text (n = 4) . The time and number scales are arbitrary and have been chosen 
appropriately in order to compare the effective simulation time ranges. The initial 
numbers of particles are respectively N = 2.25 x 10®, N = 20000 and N ~ 1500 < 
1800 due to mergings right after a renormalization step [6] [9] (respectively 4, 10 and 
10 samples), (a) fits to N(t) = Nq(1 +t/to)~^ and the effective values for £. (b) 
tN(t) is now plotted v s log(l-M). For clarity, the curve for the Kirchoff dynamics has 
been enlarged by a factor 3, due to the much smaller time range, and the different 
curves have been shifted. 

Fig. 2: Normalized vortex radii distributions obtained from the diffusing vortex simulation 
for different times. This is compared to the exact mean-field one (see text). 

Fig. 3: For the Batchelor case (n = 2), N(t) is shown for the diffusing vortex model (40 
samples) and Kirchoff dynamics (10 samples), (a) In the same units, N(t) is plotted 
vs (l+t). 
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